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FAST  ION  STOPPING  POWER  IN  DENSE,  IONIZED  PLASMAS 


I.  INTRODUCTION 

The  growing  interest  in  ion  beams  as  drivers  in  Inertial  Confinement 
Fusion  (ICF)  has  created  a  need  for  reliable  methods  of  calculating  beam 
stopping  power  as  a  fimction  of  the  temperature,  density,  and  composition  of 
the  target  material.  Such  Information  Is  essential  for  the  optimal  design  of 
ICF  target  pellet  configurations.  Because  the  slowing  down  process  can 
generate  energetic  x-rays,  it  is  also  of  interest  In  simulation  studies  for 
nuclear  weapons  effect.  It  is  therefore  evident  that  ion  stopping  power 
studies  are  useful  In  a  number  of  areas  of  Interest  to  the  scientific  and  the 
military  communities. 

In  this  report ,  a  "local  oscillator"  model  for  bound  electron  stopping 
power  is  described  and  applied  with  good  success  to  a  number  of  cold  target 
calculations.  This  model  is  then  extended  to  partially  ionized  atoms  and  com 
bined  with  free  electron  stopping  power  equations  to  calculate  total  electron 
stopping  power  results  for  target  plasmas  of  given  density,  temperature,  and 
degree  of  ionization.  Results  are  compared  with  other  calculations. 

II.  BOUND  ELECTRON  STOPPING  POWER 

Two  widely  used  models  in  the  calculation  of  the  stopping  power  of  bound 
atomic  electrons  are  the  LSS  model^,  which  is  valid  at  low  velocities,  and 
the  Bethe  theory,  which  applies  to  high  velocity  projectiles.  These  are  com¬ 
plementary  models  and  can  be  used  to  span  the  entire  velocity  range;  for 
2 

example,  Mehlhom  computes  results  for  both  models  at  a  given  velocity  and 

takes  the  lower  of  the  two  values  as  the  stopping  cross  section. 

Maniucript  submitted  January  28, 1981. 


The  chief  difficulty  In  applying  the  Bethe  theory  Is  obtaining  the 
average  Ionization  potential  I.  A  calculation  of  this  quantity  Is  a  tedious 


exercise;  hence,  various  scaling  mechanisms  have  been  proposed  to  estimate 
—  r  3  2I 

I  e.g..  Mosher  ,  Mehlhom  1  ;  to  circumvent  this  difficulty,  the  "local 

r 

oscillator  model"  (LOM)  I  Nesbet  and  Ziegler  J  was  chosen  to  represent  stopping 
by  bound  electrons.  In  the  LOM,  the  bound  electron  stopping  cross  section  Is 
given  by 


S^(v) 


4it  e^ 
mV^ 


/  4,  r= 


D(r)  K^(t)  dr. 


(1) 


where  is  the  effective  charge  of  the  projectile  ion,  e  and  m  are  the  elec¬ 
tron  charge  and  mass,  V  Is  the  projectile  velocity,  c(r)  is  the  local  bound 
electron  density  in  the  target  atom,  and  is  a  modified  zero  order  Bessel 

function. 

T  •  -h  ui  (r)/mV^,  (2) 

o 

where  .<1  is  Planck's  constant  divided  by  2ir ,  and  u)^(r)  is  the  local  plasma 
frequency  at  radius  r  within  the  atom. 

4 

There  are  three  assumptions  underlying  the  LOM  .  The  first  is  that  a 
loss  function  can  be  defined  dependent  only  on  the  local  electron  density  in 
the  target  atom.  The  second  is  that  the  longitudinal  dielectric  response  can 
be  represented  by  e(u))  with  a  single  zero  at  u  subject  to  the  high 

frequency  condition 

£(oj)  1  -  w  ^/(4)^  (3) 

o 

appropriate  to  free  electrons.  The  third  assumption  is  that  the  induced 
polarization  charge  is  spread  out  from  the  ion  trajectory  to  some  finite  radius 
of  the  order  of  the  de  Broglie  wavelength  fi/mV.  This  last  assumption  is 


Justified  by  the  adiabatic  argument  that  energy  will  be  transferred  only  to 
electrons  with  ^^eloclty  less  than  V.  A  wave  packet  with  momentum  of  order 
mV  would  have  a  spatial  spread  of  order  fi/mV. 

The  electron  density  p(r)  Is  taken  from  the  Thomas-Ferml  (TF)  model 
proposed  by  Zlnk^.  For  an  Isolated,  zero-temperature  atom,  the  TF  potential 
Is  assumed  to  have  the  form 

-  Ze/r(l  +  Ar)^,  (4) 

where  the  parameter  A  Is  determined  by  the  condition 

0^(r)  4iTr^  dr  »  Z,  (5) 

o 

p^(r)  Is  given  by  the  TF  equation 

'!<'>■  (^) 

and  Z  Is  the  number  of  electrons  bound  to  the  nucleus.  The  solution  Is 

A  -  1.14  X  10®  Z^'"®  cm‘^.  (7) 

This  leads  to  the  following  expression  for  the  TF  electron  density  of  an 
Isolated  atom  at  zero  temperature; 

c(r)  -  -^ZA^/Zirr  (1  +  Ar)^.  (8) 

As  shown  by  Zink®,  this  procedure  can  be  applied  to  a  partially  Ionized 
atom.  The  atomic  volume  Is  divided  Into  two  regions.  The  Inner  region 
around  the  nucleus,  bounded  by  r^,  contains  the  bound  electrons  and  has  a 
high  electron  density;  the  outer  region  contains  free  electrons  and  has 

a  low  density.  The  outer  region  Is  defined  by  the  Interval  where 

R  Is  the  radius  of  an  unionized  atom  at  the  given  density  of  the  material.  The 
electron  density  Is  given  by 


P  (r) 


*  2 
_Z _  .  3Z 

4  3^  4  * 

y ir  R  2irr(l+AR) 


0  £  r  ^ 


P  (r) 


£  r  £  R 


(9) 


where  Z  Is  the  number  of  free  electrons. 

Zink^  establishes  the  following  relationship  between  Z*  and 


Z* 


2  Z  Ar, 


(l+Ar^^)- 


(1+Arj^)‘ 


(10) 


Thus,  if  either  Z*  or  r^^  is  known,  the  TF  electron  density  can  be  written 
for  the  ionized  atom. 

The  approach  taken  in  this  calculation  is  to  solve  an  LIE  Saha  model  for 
the  average  charge  Z  of  the  atoms  as  a  function  of  density  and  electron 
temperature.  Then  Z*  is  set  equal  to  Z,  and  a  TF  electron  density  can  be 
written  for  an  "average  atom"  in  the  plasma.  In  this  way,  the  bound  electron 
stopping  power  in  a  partially  ionized  plasma  can  be  computed. 

The  relation  given  by  Brown  and  Moak^  is  used  to  define  the  effective 
charge  Z^  of  the  projectile  atom  with  nuclear  charge  Zp; 


Zj^/Zp  -  1,0  -  1,034  exp  (-V/V^), 


(11) 


where 

-  Zp°'^^  (e^M),  (12) 

This  relation  was  obtained  by  plotting  several  hundred  data  points  for  effec¬ 
tive  charge  and  performing  a  least-squares  fit  to  the  results.  It  was 

found  empirically  in  applying  the  LOH  model  to  cold  targets  that  good  agree- 

7  8 

ment  with  standard  results  ’  was  obtained  with  Z^  replaced  by  Z  ,  the  pro- 

0 

Jectile  atomic  number  in  Eq.  (1),  and  the  result  scaled  by 


-4- 


f  - 


(13) 


(■.h/  , 

where  y  Is  the  result  of  applying  Eq.  (11)  to  the  Incident  Ion,  and  y  Is  the 

P 

2 

result  of  Eq.  (11)  for  a  proton  at  the  same  velocity  (with  »  e  /fi) .  For 
Incident  protons,  Zeff/Z  *  1.0  gave  good  results  for  the  energies  considered. 

The  Integral  In  Eq.  (1)  Is  taken  over  the  radius  of  the  atom  (R  for 
neutral  atoms  and  r^  for  Ionized  atoms).  Low  energy  projectiles  do  not  pene¬ 
trate  the  atomic  electron  cloud  very  deeply;  hence,  only  a  few  points  near  the 
outer  edge  of  the  radial  grid  contribute  significantly  to  the  Integral.  Also, 
In  reality,  the  radius,  or  boundary,  of  the  atom  Is  not  sharply  defined.  At 
low  velocity,  from  Eq.  (11),  the  projectile  atom  has  a  small  effective  charge; 
It  retains  most  of  Its  own  electron  charge  cloud  which  Increases  the  size  of 
the  Interaction  region  with  the  target  atom.  It  Is  therefore  possible  to 
omit  significant  contributions  to  ths  stopping  cross  section  at  low  velocities 
by  using  Eq.  (1).  Indee.l,  It  sometimes  occurs  chat,  at  low  energy,  the  results 
given  by  Eq.  (1)  are  coo  low  as  compared  with  accepted  values  for  stopping 
cross  sections  e.g.,  Andersen  and  Ziegler^  .  Thus,  there  may  be  a  physical 

and  numerical  difficulty  with  Eq.  (1). 

9 

According  to  Land  and  Brennan  .  for  velocities  below  the  stopping 
cross  section  car.  be  taken  to  be  linear  In  velocity.  Therefore.  In  this 
regime,  a  scaled  cross  section  Is  calculated  from 

V  <  V^.  (14) 

where  (y^)  Is  calculated  from  Eq.  (1).  The  stopping  cross  section  Is 
Chen  taken  to  be  the  larger  of  the  results  of  Eq .  (1)  or  Eq.  (14). 

For  velocities  larger  than  the  projectile  has  a  higher  effective 
charge  and  penetrates  more  deeply  Into  the  target  atomic  electron  cloud, 
and  the  smeared  out  boundary  problem  does  not  arise. 
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Some  results  from  this  model  for  stopping  cross  sections  for  hydrogen 

atoms  In  cold,  solid  density  targets  are  shovm  In  Figure  1  and  compared  with 

the  results  of  Andersen  and  Ziegler^.  In  Figure  la,  the  agreement  Is  very 

good  over  the  entire  energy  range  for  an  kl  target.  In  Figure  lb,  for  an  Fe 

target,  the  agreement  Is  good,  except  In  the  100-200  keV  range,  where  the  LOM 

result  Is  about  20  percent  low.  For  the  Au  target  In  Figure  Ic  the  errors  are 

about  14  percent  or  less.  The  LOM  peaks  appear  to  be  shifted  to  lower  energy 

by  about  20  keV  to  50  keV  relative  to  the  reference  curves.  In  the  high 

energy,  or  Bethe  theory,  regime,  the  agreement  is  very  good.  As  noted  earlier, 

in  the  LOM  model,  a  calculation  of  mean  Ionization  potential  Is  not  necessary. 

In  Figures  Id  to  If  are  shown  comparisons  of  range  calculations  vs. 

energy  from  this  model  with  Andersen  and  Ziegler^  results.  Good  agreement  is 

obtained  again  except  at  the  lowest  energies;  this  results  from  an  omission  of 

energies  below  1  keV  In  these  calculations. 

Cold  target  calculations  for  C  atom  projectiles  are  shown  in  Figure  2. 

In  Figures  2a  and  2b,  comparisons  with  the  results  of  Northcllffe  and 
8 

Schilling  for  ki-  and  Au  targets  at  solid  density  are  given;  good  agreement  is 
obtained.  Not  shown  are  the  range  calculations  for  these  situations;  they 

g 

also  agree  well  with  Northcllffe  and  Schilling  .  Coincidences  of  the  peaks 
may  be  noted  in  Figure  2. 

The  previous  cases  dealt  with  light  ion  projectiles  on  heavier  atom 

targets.  The  situation  treated  in  Figures  3a  and  3b  is  the  reverse;  here 

cold  target  calculations  at  solid  density  are  presented  for  U  atoms  impinging 

on  A;  and  Au  atoms.  For  comparison,  the  results  of  Northcllffe  and 
8  10 

Schilling  and  Brueckner  and  Metzler  are  also  given.  In  Figure  3a,  at  low 

g 

energy,  the  LOM  results  agree  very  well  with  Northcllffe  and  Schilling  but 
are  much  lower  than  Brueckner  and  Metzler^^.  As  energy  increases,  however. 
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the  Northcliffe  and  Schilling  results  drop  much  faster  than  the  LOM  results, 
while  the  Brueckner  and  Metzler^*^  curve  approaches  the  LOM  results.  In 
Figure  3b,  the  Bruckner-Metzler^^  results  and  the  LOM  agree  very  well  at  all 

g 

energies,  while  the  Northcliffe  and  Schilling  curve  is  much  lower.  As  argued 

10  8 

by  Brueckner  and  Metzler  ,  the  Northclif fe-Schilling  data  for  heavy  ions  at 
high  energy  are  based  on  extrapolations  from  low  energy  data  and  are  therefore 
somewhat  uncertain.  Brueckner  and  Metzler^*^  calculate  the  bound  electron 
stopping  power  by  using  a  Thomas  Fermi  model  to  evaluate  the  electronic  exci¬ 
tation  energies  and  the  Coulomb  logarithm  term  in  the  stopping  power  equation; 
hence,  there  is  some  similarity  between  their  model  and  the  LOM.  However, 
they  use  a  different  effective  charge  scaling  law  which  results  in  a  higher 
effective  charge. 


III.  FREE  ELECTRON  STOPPING  POWER;  PLASMA  TARGETS 

For  heated  target  materials,  the  atoms  become  ionized,  and  stopping  due 
to  free  electrons  must  be  considered.  The  free  electron  stopping  cross  section 


is  calculated  from 


S^(V)  = 


27re^  Z 


+  An  (1  +  4  X  ) 


where  Z  is  the  average  charge  of  the  target  ions  (currently  obtained  from  a 
Saha  LTE  equation-of-state) ,  D  is  the  Debye  shielding  length,  and 

X  -  (mV^/2kT)^^^,  (16 


where  kT  is  the  electron  temperature. 


F(x)  =  erf(x)  -  2  x  e  ^ 
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The  quantity  b  is  the  minimum  impact  parameter  for  electron-ion 


scattering  and  is  given  by 

i 

b  -  MAX  ' 


mV 


(18) 


i.e.,  the  maximum  of  either  the  classical  or  quantum-mechanical  impact  para¬ 
meter  defined  by  the  uncertainty  principle. 

The  first  term  in  Eq.  (15)  is  the  short. range  ion-electron  binary- 

encounter  scattering  term  and  is  taken  from  the  work  of  Campbell^^.  who 

12 

adapted  it  from  Brueckner  and  Brysk  .  The  second  term  is  the  polarization 

13 

term  and  is  taken  from  Pines  and  Bohm  .  For  distances  larger  than  D,  the 
plasma  acts  as  a  continuous  medium,  and  distant  collisions  cause  loss  of 
energy  by  the  excitation  of  plasma  oscillations,  which  appear  as  an  oscil¬ 
lating  wake  behind  the  projectile. 

2 

.As  discussed  by  Mehlhom  ,  the  effective  charge  of  the  projectile  ion 
will  be  increased  by  the  plasma  free  electrons  due  to  increased  collisional 
ionization.  A  higher  relative  velocity  occurs  between  ion  and  plasma  than  for 
the  cold  target.  Therefore,  in  Eqs.  (11)  and  (18),  the  ion  velocity  is  re¬ 
placed  by  the  addition  of  the  ion  velocity  and  electron  thermal  velocity  in 
a  random  phase  manner.  The  scaling  defined  by  Eq.  (13)  is  also  applied  to 
Eq.  (15). 

In  a  plasma  target,  the  bound  atomic  electrons  can  be  screened  from  the 

projectile  ions  by  the  free  electrons.  Thus,  the  LOM  must  be  modified  to  take 

this  shielding  effect  into  account.  For  an  electron  in  an  isolated  atom,  the 

maximum  impact  parameter  is  V/u,  where  oj  is  a  characteristic  frequency  of 
14  x5 

motion  ’  .  Plasma  screening  limits  this  parameter  to  the  Debye  length  D. 

Eq.  (2)  can  be  rewritten 
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iA. 


— - ^ 


T  »  (h/mV)  *  (o)  (r)/V)  . 

o 

Whenever  V/u  (r)  >  D,  t  is  taken  to  be 
o 

T  =  h/mVD  (19) 

Thus,  plasma  shielding  of  the  bound  electrons  is  taken  into  account  by  limiting 

the  argument  of  the  Bessel  function  as  given  by  Eq.  (19). 

For  a  plasma  target  at  a  given  temperature  and  ion  density,  a  Saha  LTE 

model  is  used  to  obtain  the  average  ionic  charge  Z.  Using  the  LOM  as  modified 

for  plasma  screening  in  combination  with  Zink's^  model  for  determining  Thomas- 

Fermi  electron  densities  for  partially-ionized  atoms  gives  the  bound  electron 

stopping  cross  section.  Adding  this  result  to  Eq.  (15)  gives  the  total 

electronic  stopping  cross  section. 

This  formalism  has  been  applied  to  the  stopping  power  of  Au  at 

kT  =  200  eV  and  at  0.01  times  solid  density  for  carbon  ions.  The  results 

2 

are  shown  in  Figures  4a  and  4b;  also  given  are  Mehlhom's  results  for  com¬ 
parison.  Mehlhom's  S,  (V)  results  are  generally  higher,  especially  at  high 

0 

energy,  where  the  difference  is  'v  25  percent.  As  mentioned  earlier, 

Mehlhom  uses  a  Bethe  model  to  calculate  bound  electron  stopping  at  high 
energy;  the  differences  could  be  due  to  different  effective  I's  in  the  two 
calculations.  Also,  Mehlhom  ignores  plasma  shielding  effects  on  the  bound 
electrons,  which  would  lead  to  smaller  values  of  S^(V).  Figure  4b  shows 
the  free  electron  stopping  cross  sections;  here  the  agreement  is  good  over 
the  entire  energy  range.  Since  free  electron  stopping  clearly  dominates  the 
bound  electron  contribution  in  this  case,  the  total  electronic  stopping  power 
is  approximately  the  same  for  both  calculations,  especially  at  the  higher 
energies . 

Figure  5  shows  a  comparison  with  Brueckner  and  Metzler^^  for  H  atoms  in 
an  ki  target  at  solid  density  and  at  kT  =  10  eV  and  2  keV.  At  2  keV,  A£  is 
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completely  Ionized,  and  the  electronic  stopping  Is  totally  due  to  free  elec¬ 
trons.  The  two  calculations  are  In  good  agreement  here.  At  10  eV,  bound 
electron  stopping  cross  sections  are  larger  than  the  free  electron  cross 
sections  by  about  50  percent;  hence  screening  of  the  bound  electrons  Is 
Important  In  this  case.  The  Brueckner-Metzler  result  Is  about  30  percent 
lower  than  the  results  of  this  model.  As  mentioned  earlier,  they  use  a 
different  model  to  calculate  bound  electron  stopping.  They  take  plasma 
screening  Into  account  by  using  a  combination  of  the  Debye  length  D  and  the 
Isolated  atom  maximum  Impact  parameter  V/u. 

Figure  6  shows  total  electronic  stopping  cross  sections  for  protons  on 
solid  density  Au  at  kT  =  10  eV  and  2  keV;  also  shown  are  the  Brueckner- 
Metzler  calculations.  At  10  eV,  stopping  by  bound  electrons  dominates  the 
cross  sections;  the  two  calculations  are  In  very  good  agreement.  At  2  keV, 
free  electron  stopping  predominates.  The  Brueckner-Metzler  results  are  about 
20  to  25  percent  lower;  this  Is  probably  due  to  differences  In  the  models  for 
free  electron  stopping  power. 

Figure  ’  shows  cross  section  calculations  for  U  Ions  on  a  solid  density 
All  target.  At  kT  =  10  eV,  there  Is  a  wide  disparity  between  these  results 
and  the  Brueckner-Metzler  calculations  at  low  energy;  at  higher  energy,  their 
results  are  approximately  30  percent  lower.  At  2  keV,  there  Is  a  large  dis¬ 
agreement  In  the  low  energy  range,  but  the  results  appear  to  be  merging  at 
the  high  energy  end.  These  differences  appear  to  be  explained  by  the  differ¬ 
ence  In  the  effective  charge  scaling  laws.  According  to  Brueckner  and  Metzler, 

2 

their  scaling  law  can  show  a  ratio  (Z^^^/Z)  as  large  as  twice  the  result  of 
Eq.  (11)  for  energies  of  few  MeV/amu. 

Figure  8  shows  comparisons  with  Brueckner  and  Metzler  for  U  Ions  at  one- 
tenth  solid  density  Au  at  kT  “  200  eV  and  2  keV.  The  2  keV  results  are 
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similar  In  appearance  to  those  of  Figure  7,  and  the  same  statements  regarding 
scaling  laws  apply.  The  200  eV  curves  are  In  reasonably  good  agreement  for 
energies  ^  4  MeV/amu.  Free  electron  stopping  dominates  by  factors  of  1.5  to 
2  at  higher  energies. 

Nardi,  Peleg,  and  Zinamon^^  (NPZ)  have  calculated  electronic  stopping 
cross  sections  for  protons  in  Au;  comparisons  with  these  results  are  shown  in 
the  next  two  figures.  Figure  9  shows  free  electron  stopping  cross  sections  at 
solid  density  and  kT  =  1  keV.  The  NPZ  results  for  both  collisional  and  non- 
colllsional  plasma  dielectric  functions  are  shorn;  also  given  are  calculations 
based  on  Eq.  (15)  and  the  Brueckner-Metzler  results.  The  agreement  between 
this  model  (Eq.  (15))  and  the  NPZ  non-collisional  results  is  good  over  the 
energy  range  1  to  10  MeV.  The  Brueckner-Metzler  curve  agrees  with  the  NPZ 
collisional  results  at  low  energy  but  falls  below  the  other  calculations  as 
energy  increases;  at  10  MeV,  the  Brueckner-Metzler  result  is  approximately 
20  percent  lower  than  the  others. 

Eq.  (15)  represents  a  type  of  model  that  NPZ  refer  to  as  "binary  plus 
collective",  i.e.,  short  range  binary  collisions  plus  long  range  crllective 
effects;  at  high  energies  (^  4  MeV),  they  obtained  good  agreement  between  this 
kind  of  model  and  their  more  complex  dielectric  function  models.  The  results 
of  Eq.  (15)  as  shown  in  Fig.  9  reinforce  their  conclusions;  in  fact,  the 
Eq.  (15)  results  are  in  excellent  agreement  with  their  "binary  +  Pines-Bohm" 
result.  At  low  energies,  where  the  projectile  velocity  is  not  large  compared 
to  the  electron  thermal  velocity,  the  dielectric  function  models  predict  lower 
cross  sections  than  the  binary  plus  collective  models. 

Figure  10  shows  free  electron  stopping  cross  sections  for  protons  in  Au 

at  one-hundredth  of  solid  density  and  kT  ■  1  keV.  The  NPZ  dielectric  function 

2 

results,  calculations  based  on  Eq.  (15),  and  Mehlhom's  dielectric  function 

14 

and  binary-plus-collective  (i.e.,  Jackson  )  results  are  given.  The 
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dielectric  function  calculations  of  NPZ  and  Mehlhom  agree  -ery  well  in  the 
1  to  10  MeV  range.  Mehlhom 's  blnary-plus-collectlve  results  lie  between  the 
dielectric  function  calculations  and  the  results  from  Eq.  (15).  Again  there 
is  some  disparity  at  lower  energies,  but  all  the  curves  tend  to  merge  at  high 
energy. 

Figure  11  shows  free  electron  stopping  cross  sections  for  protons  in  Au 

at  one-tenth  solid  density  and  kT  >»  50  eV.  Results  from  Eq.  (15)  are  shown 

2 

along  with  Mehlhom 's  dielectric  function  calculations  and  his  binary  plus- 
collective  calculations.  Agreement  is  good  among  the  three  results  for  energy 
^  3  MeV.  This  again  bears  out  the  conclusion  of  NPZ  and  also  Mehlhom  that, 
in  this  energy  range,  the  binary-plus-collective-oscillations  model  is  in 
good  agreement  with  more  complex  dielectric  function  calculations. 

The  dependence  of  this  model  on  electron  temperature  is  displayed  in 
Figures  12  and  13.  In  both  calculations,  the  cross  sections  at  10  eV  are 
lower  than  the  cold  target  cross  sections;  this  is  due  to  shielding  of  the 
bound  electrons  by  the  free  plasma  electrons.  As  the  electron  temperature 
rises  and  the  free  electron  contribution  becomes  more  significant,  the  cross 
section  rises  again  so  that,  at  50  eV,  they  are  larger  than  the  cold  target 
cross  sections.  At  energies  above  5  MeV/amu  for  these  calculations;  the 
cross  sections  continue  to  rise  with  increasing  kT.  At  lower  energies  as 
kT  rises,  the  cross  sections  show  a  decrease.  This  is  due  to  the  function 
r'(x)  (defined  by  Eq.  (17)).  This  function  is  a  measure  of  the  Coulomb  inter¬ 
action  between  the  projectile  ion  and  the  plasma  electrons  and  restricts  energy 
losses  by  the  projectile  ion  only  to  electrons  with  velocity  lower  than  the 
ion  velocity.  As  kT  Increases,  losses  by  low  velocity  ions  to  the  plasma 
decrease. 

This  initial  drop  in  the  stopping  power  as  the  temperature  rises  in  the 
target  material  is  the  same  trend  predicted  by  Brueckner  and  Metzler^^; 
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however,  their  model  predicts  a  much  more  severe  decrease  than  this  model 

indicates.  This  is  probably  due  to  differences  in  the  models  and  in  different 

treatments  of  the  plasma  shielding  effects. 

Calculations  with  this  model  which  omit  this  screening  effect  predict  a 

2 

steady  increase  of  stopping  power  with  electron  temperature,  as  Mehlhom  pre¬ 
dicts;  there  is  no  initial  drop  as  kT  rises  above  the  cold  target  temperature  . 
Mehlhom  omits  this  plasma  screening  effect. 

It  should  be  noted  that  the  differences  between  the  cold  target  and  the 
1  keV  cross  sections  in  Fig.  12  are  in  good  agreement  with  the  results  given 
by  NPZ^^  for  their  cold  target  and  "binary  +  Pines  Bohm". 

The  effect  of  density  in  this  model  is  shown  in  Fig.  14,  where  stopping 

cross  sections  for  U  ions  in  Au  and  solid  density  and  0.1  p^  are  given  for 

10  eV,  100  eV,  and  1000  eV.  In  each  case,  the  lower  density  target  gave 

higher  stopping  cross  sections.  This  effect  is  also  predicted  by  Brueckner 

10  2 
and  Metzler  and  Mehlhom  . 

3 

Mosher  has  calculated  a  dimensionless  correction  term  (In  A  -  in  A  ) 

s 

which  is  a  measure  of  the  difference  between  stopping  powers  for  a  heated 
plasma  and  for  a  cold,  neutral  target;  he  shows  the  variation  of  this  quantity 
for  targets  with  Z  ■  6  and  Z  *  81  for  the  range  10  eV  ^  kT  ^  1000  eV  at 
several  electron  densities.  In  an  effort  to  compare  with  these  results,  some 
of  the  stopping  cross  sections  calculated  here  were  taken  at  a  given  energy, 
the  cold  target  results  at  the  same  energy  were  subtracted,  and  the  difference 
divided  by  the  quantity 

Z  (20) 

3 

to  obtain  a  normalized,  unitless  difference  analogous  to  Mosher's  . 


4ire^ 

m 

e 
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were  chosen  for 


Results  for  U  projectiles  on  Au  targets  at  p  ■  0.1 
j  comparison;  for  this  situation,  Z  •  79,  and  the  product  of  Z  times  atomic 

number  density  is  4.6  x  10“  ;  this  is  reasonably  close  to  his  Z  ■  81, 
i  S  ■  10^^  case.  Calculations  were  done  for  energies  of  12  MeV/amu  and 

20  MeV/amu  for  this  model,  and  similar  calculations  were  done  for  these 
energies  from  the  Brueckner-Metzler^^  data.  The  results  are  summarized  in 
Table  1. 

TABLE  1:  U-*Au,  p-0.1c,N  -  10“^ 

_  s  e 


1  1 

_ 1 

E  ■  12  meV/amu 

E  ■  20  MeV/amu  j 

!  ! 

IkTCeV)  i 

i  1 

;  1 

Mosher^ 

This  model 

B-M^O 

1 

This  model  ' 

b-m'°  : 

!  10 

[  i 

1 

0  i 

-  0.14  i 

i 

-  0.64 

-  0.18  1 

1 

1 

-  0.69 

+  0.21 

1  -  0.37 

1 

+  0.23 

1  -  0.48 

+  1.68 

;  +  1.80 

-J _ 

1  +  1.68 

1  +  1.75 

( 

j 

3 

Mosher's  correction  term  is  independent  of  projectile  velocity  and 

depends  only  on  the  target  plasma.  Both  the  Brueckner-Metzler^^  results 

and  calculations  from  this  model  appear  to  verify  this  conclusion  in  this 

energy  range.  The  model  used  in  this  report  agrees  fairly  well  with 
3 

Mosher  at  10  eV  and  at  100  eV  and  with  Brueckner  and  Metzler  at  1000  eV, 

3 

where  Mosher's  result  is  larger  by  a  factor  of  order  2.5.  The  Brueckner- 

10  3 

Metzler  results  are  much  lower  than  Mosher's  at  each  value  of  kT. 

3 

Mosher  also  presents  results  for  this  correction  term  for  a  target  Z 

of  6  and  several  electron  densities.  Calculations  were  done  for  protons 

incident  on  carbon  atoms  at  E  ~  10  MeV/amu  and  20  MeV/amu  and  at  electron 
20  23 

densities  of  10“  and  10  .  Comparisons  are  given  in  Table  2. 
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TABLE  2:  H  -C 


1 

i 

o 

O 

10^^ 

'kT  (eV) 

Mosher^ 

amu 

E-20^ 

amu 

1  ~ 
1  Mosher’’  [ 

E-io5^£i^ 

amu 

MeV 
E-20  — 
amu 

'  10 

.0 

+  3.04 

+  3.08 

1  i 

!  1 

:  -  0.5  j 

i  1 

+  0.03 

+  0.16 

;  30 

+0.2 

+  3.25 

+  3.25 

1  1 

!  -  1.0  ! 

I  1 

+  1.1-3 

I 

+  1.06 

.  100 

+  2.8 

+  5.44 

■ 

+  5.50 

,  ^  +  0.4  ; 

+  1.86 

+  1.84 

300 

+4.5 

+  5.50 

t 

+  5.50 

f  t 

i  ”-  +  1.7  1 

J 

+  2.40 

+  2.40 

i 

1000 

+  6.4 

+5.34  i 

I 

+  5.53 

i  ^  +  2.8  1 

'  » 

_ L 

+  2.34 

+  2.44  : 

i 

The  results  of  this  model  are  again  approximately  Independent  of  velocity 

3  '*0 

In  this  energy  range.  Mosher's  results  for  -  10  vary  over  a  much  wider 

range  as  kT  changes;  this  may  be  due  to  a  difference  of  Ionization  models. 

23  3 

The  results  at  ■  10  are  generally  higher  than  Mosher's  except  at  1  keV, 
where  they  are  comparable. 

3 

The  differences  with  Mosher's  results  are  probably  due  to  different  free 
electron  stopping  power  models. 


IV.  CONCLUSIONS 

A  local  oscillator  model  for  bound  atomic  electron  stopping  power  has 

been  utilized  to  calculate  stopping  cross  sections  for  several  projectile  ions 

in  cold  targets;  Thomas  Fermi  electron  densities  as  proposed  by  21nk^  were 

7  8 

us^d  in  the  computations.  Good  agreement  with  standard  results  '  were 
obtained  over  the  entire  energy  range.  Tnus  the  LOM  appears  to  have  a 
solid  basis  as  a  model  for  bound  electron  stopping  power. 

Using  Zink's^  formulation,  this  model  was  extended  to  partially  ionized 


atoms  and  modified  to  take  into  account  shielding  of  the  bound  atomic  electrons 


by  the  free  plasma  electrons.  It  was  then  combined  with  a  binary-collisions- 
plus-collective-oscillations  model  of  free  electron  stopping  power  to  calculate 
total  electron  stopping  power  in  heated  plasmas. 

In  most  instances,  reasonable  agreement  with  comparison  calculations  was 
obtained;  hence  some  confidence  in  these  results  are  justified. 

From  an  examination  of  Figures  4  through  11,  it  is  apparent  that,  in  the 
treatment  of  free  electron  stopping  power,  this  model  is  closer  to  those  of 
Mehlhorn^  and  NPZ^^  than  it  is  to  that  of  Brueckner  and  Metzler^^.  Both  NPZ 
and  Mehlhorn  agree  chat,  in  the  high  energy  range,  the  so-called  binary-plus- 
collective  oscillations  approach  as  represented  by  Eq.  (15)  yields  results  in 
good  agreement  with  more  complex  dielectric  function  calculations.  There  are 
differences  in  the  treatment  of  bound  electron  stopping  power:  Mehlhorn  and 
NPZ  use  the  Bethe  theory,  while  this  report  adopts  the  local  oscillator  model. 

One  of  Che  questions  still  to  be  settled  is  the  effective  charge  on  the 
projectile  ion.  The  expression  obtained  by  Broxm  and  Moak^  was  used  here,  but 
Brueckner  and  Metzler^*^  adopt  an  expression  that  yields  a  higher  effective 
charge . 

.■Another  problem  which  this  model  does  not  address,  but  which  is  under 
study,  is  the  problem  of  electron  degeneracy  which  arises  at  densities  well 
above  solid  densities.  .-Iny  detailed  study  of  ICF  plasmas  must  Include 
degeneracy  effects. 
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Fig.  la  —  Cold  target  stopping  cross  section  of  solid  density  Ai;  for  H  ions 
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"^.g.  lb  —  Cold  target  shopping  cross  section  of  solid  density  Fe  for  H  ions 
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Fig.  Id  ~  Range  vs.  energy  of  H  ions  in  cold  target,  solid  density  AS. 
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Fig.  2a  —  Cold  target  stopping  cross  section  of  solid  density  A£  for  C  ions 
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Fig.  2b  —  Cold  target  stopping  cross  section  of  solid  density  Au  for  C  ions 
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Fig.  5  —  Total  stopping  cross  section  of  solid  density  A£  at 
kT  =  10  eV  and  kT  =  2  ke\’  for  H  ions 
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Fig.  6  —  Total  stopping  cross  section  of  solid  density  Au  at 
kT  =  10  eV  and  kT  =  2  keV  for  H  ions 
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Fig.  9  —  Free  electron  stopping  cross  sections  for  protrons  in 
solid  density  Au  at  kT  =  1  keV 
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Fig.  10  —  Free  elect!  ■  topping  cross  sections  for  protrons  in  Au  at 
0.01  tunes  solid  density  and  kT  *  1  keV 
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Fig.  11  —  Free  electron  stopping  cross  sections  for  protrons  in  Au  at 
0.1  times  solid  density  and  kT  •  50  eV 
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Fig.  12  —  Total  stopping  cross  sections  for  protrons  in  solid 
density  Au  at  several  values  of  kT 
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Fig.  13  —  Total  stopping  cross  sections  for  U  ions  in  Au  at  0.1  times 
solid  density  and  at  several  values  of  kT 
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